Lecture 9 - Modeling, Simulation, and 

Systems Engineering 


• Development steps 

• Model-based control engineering 

• Modeling and simulation 

• Systems platform: hardware, systems software. 
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Control Engineering Technology 

• Science 

- abstraction 

- concepts 

- simplified models 

• Engineering 

- building new things 

- constrained resources: time, money, 

• Technology 

- repeatable processes 

• Control platform technology 

• Control engineering technology 
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Controls development cycle 

• Analysis and modeling 

- Control algorithm design using a simplified model 

- System trade study - defines overall system design 

• Simulation 

- Detailed model: physics, or empirical, or data driven 

- Design validation using detailed performance model 

• System development 

- Control application software 

- Real-time software platform 

- Hardware platform 

• Validation and verification 

- Performance against initial specs 

- Software verification 

- Certification/commissioning 
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Algorithms/Analysis 

Much more than real-time control feedback computations 

• modeling 

• identification 

• tuning 

• optimization 

• feedforward 

• feedback 

• estimation and navigation 

• user interface 

• diagnostics and system self-test 

• system level logic, mode change 
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System and software Controls analysis 


Model-based Control Development 


Conceptual 

Analysis 


Control design model: 
x(t+l ) = x(f) + u(t) 




Conceptual control 
algorithm: 
u = -k(x-x d ) 


Application 
code: Simulink 


Validation and 
verification 


Deployment 


Detailed simulation 
model 


f 


< 


Hardware-in-the- 
loop sim 


c 


Physical plant 


> 


W 


Detailed control application: 

saturation, initialization, BIT, 
fault recovery, bumpless transfer 



Prototype 
controller 


Systems platform: 
Run-time code, OS 
Hardware platform 


Deployed 

controller 
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Controls Analysis 


Conceptual 

Analysis 


Data model 


Fault model 


Control design model: 
x(t+l ) = x(/) + u(t) 


<C 


o 


Identification & tuning 






Accomodation 

algorithm: 


Conceptual control 
algorithm: 
u = -k(x-x d ) 


Application 
code: 
Simulink 



Detailed 

simulation 

model 


Detailed control application: 
saturation, initialization, BIT, 
fault recovery, manual/auto 
mode, bumpless transfer, 
startup/shutdown 
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The rest of the lecture 


• Modeling and Simulation 

• Deployment Platform 

• Controls Software Development 
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Modeling in Control Engineering 


Control in a 

system 

perspective 


Control analysis 
perspective 


Measurement 

system 

Sensors 

-— 



System model 




Control 

computing 



Control 


handles 

Actuators 
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Models 


• Why spend much time talking about models? 

- Modeling and simulation could take 80% of control analysis effort. 

• Model is a mathematical representations of a system 

- Models allow simulating and analyzing the system 

- Models are never exact 

• Modeling depends on your goal 

- A single system may have many models 

- Large ‘libraries’ of standard model templates exist 

- A conceptually new model is a big deal (economics, biology) 

• Main goals of modeling in control engineering 

- conceptual analysis 

- detailed simulation 
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Modeling approaches 

• Controls analysis uses deterministic models. Randomness and 
uncertainty are usually not dominant. 

• White box models: physics described by ODE and/or PDE 

• Dynamics, Newton mechanics 

x = f(x,t ) 

• Space flight: add control inputs u and measured outputs y 

x = f(x,u,t ) 

y = g(x,u,t) 
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Orbital mechanics example 



1643-1736 


Newton’s mechanics 


- fundamental laws 


- dynamics 

r 

v — —ym — 
r 



r — v 




Laplace 

- computational dynamics 
(pencil & paper computations) 

- deterministic model-based 
prediction 


1749-1827 
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x = f(x,t ) 


X = 


Vi 


v. 


V, 







Sampled and continuous time 

• Sampled and continuous time together 

• Continuous time physical system + digital controller 

- ZOH = Zero Order Hold 


A/D, Sample 


Control 
:N computing 



D/A, ZOH 


D - 

0 


^- 



o - 







Actuators 
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Servo-system modeling 

• Mid-term problem 

• First principle model: electro-mechanical + computer sampling 

• Parameters follow from the specs 



my + j3y + b(y-x) + c(y-x) = F 
Mx + b(x — y) + c(x — y) = 0 
F = fl, Tjl +1 - gu 
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Starting Point 


Finite state 
machines 


• TCP/IP State Machine 
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Hybrid systems 


Combination of continuous-time dynamics and a state machine 
Thermostat example 

Analytical tools are not fully established yet 

^m to 

Simulation analysis tools are available 
- Stateflow by Mathworks 
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PDE models 


Include functions of spatial Example, sideways heat equation 
variables 

- electromagnetic fields 

- mass and heat transfer 

- fluid dynamics 

- structural deformations 

For ‘controls’ simulation, r 
reduction step is necessary 



- Usually done with FEM/CFD data 

- Example: fit step response 


dT , d 2 T 
— = k —=- 
dt dx 2 

T( 0) = u\ T( 1) = 0 

dT 
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Simulation 


• ODE solution 

- dynamical model: x = f(x,t ) 

- Euler integration method: x(t + d) = x(t) + d ■ f{x(t),t) 

- Runge-Kutta method: ode45 in Matlab 

• Can do simple problems by integrating ODEs 

• Issues with modeling of engineered systems: 

- stiff systems, algebraic loops 

- mixture of continuous and sampled time 

- state machines and hybrid logic (conditions) 

- systems build of many subsystems 

- large projects, many people contribute different subsystems 
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Simulation environment 



• Block libraries 

• Subsystem blocks 
developed independently \ 

• Engineered for developing., 
large simulation models 

• Controller can be designed 
in the same environment 

• Supports generation of 
run-rime control code 


Simulink by Mathworks 
Matlab functions and analysis 
Stateflow state machines 

• Ptolemeus - 
UC Berkeley 


mi 




















Model development and validation 


• Model development is a skill 

• White box models: first principles 

• Black box models: data driven 

• Gray box models: with some unknown parameters 

• Identification of model parameters - necessary step 

- Assume known model structure 

- Collect plant data: special experiment or normal operation 

- Tweak model parameters to achieve a good fit 
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First Principle Models - Aerospace 




beta 

(deg) 


P= < 

(deg/sec) 


Aircraft models 

Component and 
subsystem modeling and 
testing 

CFD analysis 

Wind tunnel tests - to 
adjust models (fugde 
factors) 

Flight tests - update » 
aerodynamic tables and 
flight dynamics models 


(deg/sec) 



+ flight k 
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time (sec) 

NASA Langley - 1998 
HARV - F/A-18 



Airbus 380: $13B development 
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Step Response Model - Process 


• Dynamical 
matrix 
control 
(DMC) 

• Industrial 
processes 


FEED 




control inputs 



Dynamic Matrix for Propvlene/Propane Splitter 


.Ell 
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Approximate Maps 

• Analytical expressions are rarely sufficient in practice 

• Models are computable off line 

- pre-compute simple approximation 

- on-line approximation 

• Models contain data identified in the experiments 

- nonlinear maps 

- interpolation or look-up tables 

- AI approximation methods 

• Neural networks 

• Fuzzy logic 

• Direct data driven models 
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Empirical Models - Maps 

• Aerospace and automotive - have most developed 
modeling approaches 


• Aerodynamic tables 

• Engine maps 

- turbines-jet engines 

- automotive - ICE 


Example 



Rolling-Moment Coefficient as a Function of 
Angle of Attack with TEF Deflections, P = 10°. 


TEF=Trailing Edge Flap 
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Empirical Models - Maps 


Process control mostly uses empirical models 


Process maps in 

semiconductor 

manufacturing 

Epitaxial growth 

(semiconductor 

process) 

- process map for 
run-to-run control 



Radius 


H 2 Flow 


"° 5 Outer Gas Flow Valve 
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Multivariable B-splines 


• Regular grid in multiple variables 

• Tensor product of B-splines 

• Used as a basis of finite-element models 


y(u,v) = ^ w j,k B j ( u )B k (v) 
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Neural Networks 


Any nonlinear approximator might be called a Neural Network 

- RBF Neural Network 



- Polynomial Neural Network 

- B-spline Neural Network 

- Wavelet Neural Network 

MPL - Multilayered Perceptron 

- Nonlinear in parameters 

- Works for many inputs 


Linear in parameters 


y(x) = w, 0 +f 


f 0) = 


f 


2N jy) Z w 2 jXj 


V j J 

l-e~ x 




V j 


J 


1 + e 


-x 


i 

y 

i —— 

/ V=f[\ 
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Layer L 
(output layer) 


Layer L-1 
(hidden layer) 



Layer 2 
(hidden layer) 


Layer 1 
(Input layer) 
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Multi-Layered Perceptions 


• Network parameter computation 

- training data set 

- parameter identification 

y(x) = F(x;9 ) 

• Noninear LS problem 

V = - F(x ij) ;0) —> min 

j 

• Iterative NLS optimization 

- Levenberg-Marquardt 

• Backpropagation 

- variation of a gradient descent 
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Neural Net application 


Internal Combustion Engine 
maps 

Experimental map: 

- data collected in a steady state 
regime for various combination 
of parameters 

- 2-D table 

NN map 

- approximation of the 
experimental map 

- MLP was used in this example 

- works better for a smooth 
surface 



Torque [Nm] 
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Fuzzy Logic 


• Function defined at nodes. Interpolation scheme 

• Fuzzyfication/de-fuzzyfication = interpolation 

• Linear interpolation in 1 -D 



very pale right too brown black 

pale colour 

• Marketing (communication) and social value 

• Computer science: emphasis on interaction with a user 

- EE - emphasis on mathematical analysis 
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Local Modeling Based on Data 


Data mining in the loop 
Honeywell Prague product 


Heat Loads 




Multidimensional 
Data Cube 



11339 J|F.b,u.„ 100:00:00 U V Actual 


1 


Data: |Steamdeivew-Aveiaqe[l/hl Model: Temperature onlv; 1 hour;dataslarlinaall.9.96 -3 Jjltl 

<-1 Show values: |121 - |168 from |1BS foiecasted ■+ | lable | 

Dale: 123 02.1999 07:00:00 Value: |30«9 Confidence: ITocT 


Heat 




demand 


Query point 
( What if ? ) 


Outdoor 
temperature 


Time 
of day 


Forecasted 

variable 


Explanatory 

variables 
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System platform for control computing 


• Workstations 

- advanced process control 

- enterprise optimizers 

- computing servers 
(QoS/admission control) 

• Specialized controllers: 

- PLC, DCS, motion controllers, 
hybrid controllers 
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System platform for control computing 








Embedded: |iP + software 
DSP 


FPGA 



ASIC / SoC 
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Embedded 


processor 

range 
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System platform, cont’d 


• Analog/mixed electric circuits 

- power controllers 

- RF circuits 


• Analog/mixed other 


EM = 

Electr-opt 

Modulator 


- Gbs optical networks 


Signal 
input 

Data {V| -► 10-Gbit/s 
EM 

Data bar (Vlh> driver IC 


Automatic 


Temperature 


temperature 


monitor 

T 

control 


interface 



ti 


EM laser: 


Optical 

output 


Laser 

temp. 

monitor 

output 




-o- 

Monitor T 
output 


T 


Shutdown 
signal (L,— 
input 


H 


L 


Shutdown 

interface 


Automatic 
► power 
control 


Alarm/ 

monitor 

interface 


XJ 


Transmitter 


Alarm/ 

monitor 

output 


pin-PD 

pre-amplifier 


AGC 

amplifier 


Optical 

input 


VCO 

jr 

Clock and 
data 

regenerator 


LOS alarm 
interface 


Data 
0 output 

V Clock 
output 

Loss of 

Y signal 
alarm 
output 


Receiver 


Functional Block Diagram of 10-Gbit/s Optical Transmitter/Receiver. AGC = Auto Gain Control 
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Control Software 


• Algorithms 

• Validation and Verification 
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Controller 
Requirements 


The Process 


Executable Algorithm Model 
Test Cases / Expected Results 



Target Processor Code 
Default Calibration 


control system 
calibration 


Calibration 

Analytical 

Dyno 

Vehicle 



New/Refined Requirements 
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Control 

application 

software 

development 

cycle 

• Matlab+toolboxes 

• Simulink 

• Stateflow 

• Real-time Workshop 
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Engineering 

Problem 


Measured Data 


Data Modeling, 
Analysis, and 
Visualization 



Behavioral 

Modeling 





Control System 
Design and 
Analysis 


Simulation 


Code Generation 


EMBEDDED 

CODE 



Hardware-in- 

the-Loop 

Simulation 



Embedded System 


DEPLOYMENT 

























Real-time Embedded Software 


Mission critical 

RT-OS with 
hard real-time 
guarantees 

C-code for each 
thread generated 
from Simulink 

Primus Epic, 
B787, A380 


Application < 


Deos 


Hardware 


Process 1 


Thread 


Thread 


Process 2 


Thread 


Thread 


Process 3 


Thread 


Thread 


A 

_ ft _ 

& v 

iL Platform 

Kernel ^ ypRegistry 

HAL 

PAL 

■v> 

CPU 

l)r 

Platform Hardware 



Deos source code 

translated for automatic analysi$\switch 
and verification of time 
partitioning property 
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Hardware-in-the-loop simulation 


Aerospace 
Process control 
Automotive 
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Requirement Specification 


Algorithm Specifications 


Virtual Prototyping 



Physical Prototyping 


Cadence 





































